Critical Entropy of Quantum Heisenberg Magnets on Simple-Cubic Lattices 
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We analyze the temperature dependence of the entropy of the spin- 1/2 Heisenberg model on the 
three-dimensional simple-cubic lattice, for both the case of antiferromagnetic and ferromagnetic 
nearest neighbor exchange interactions. Using optimized extended ensemble quantum Monte Carlo 
simulations, we extract the entropy at the critical temperature for magnetic order from a finite-size 
scaling analysis. For the antiferromagnetic case, the critical entropy density equals 0.341(5)fcs , 
whereas for the ferromagnet, a larger value of 0.401(5)fes is obtained. We compare our simulation 
results to estimates put forward recently in studies assessing means of realizing the antiferromagnetic 
Neel state in ultra-cold fermion gases in optical lattices. 



(N 



i 



i 

C 

O 

o 



> 
O 
0\ 
On 
CO 

CN 

o 



I. INTRODUCTION 

Ultra-cold atom gases provide a unique opportunity 
to study basic models of strongly interacting many-body 
systems in a highly controlled experimental framework 1 . 
After the pioneering realization of the superfluid-to-Mott 
insulator transition of ultra-cold bosonic atoms^, re- 
cently different groups reported the realization of a Mott- 
insulating state also for ultra-cold 40 K fermionic atoms on 
three-dimensional simple-cubic lattices 3, A next major 
step would be the realization and identification of a low- 
temperature antiferromagnetically ordered Neel state in 
such systems. In fact, various recent studies explored the 
prospects of realizing the antiferromagnetic Neel state in 
an ultra-cold gas of fermionic atoms confined to an op- 
tical lattice^i 2 . A quantity that is important in assess- 
ing the relevant temperature scale is the upper bound 
on the entropy allowed in order to transfer the atoms 
into the antiferromagnetically ordered state. For a three- 
dimensional fermionic Hubbard model on a simple-cubic 
lattice at half-filling, estimates of the entropy S as a func- 
tion of temperature, as well as the ratio U/t between the 
onsite repulsion U and the hopping amplitude t, were ob- 
tained within a single-site dynamical mean-field theory 
(DMFT) approach^. As reported in Ref. H, DMFT how- 
ever overestimates the Neel temperature by about 50% 
in the intermediate coupling regime, as compared to di- 
rect quantum Monte Carlo simulations based on systems 
with up to 10 3 lattice sites^ 3 .. Obtaining in a similar way 
the entropy from direct quantum Monte Carlo simula- 
tions is challenging, and typically involves integration of 
the specific heat over an extended temperature range. 

However, in the limit of large U/t, the spin physics 
of the Hubbard model is well known to be described 
by a nearest-neighbor spin-1/2 quantum Heisenberg 
model, with an antiferromagnetic exchange coupling J = 
4t 2 /U > 0, obtained within perturbation theory around 
the strong coupling limit t = 0. This model is accessi- 
ble to large scale quantum Monte Carlo simulations and 
moreover - as shown below - it is possible to provide 
precise values of the entropy in this large- U limit. In 
particular, one can obtain the value of the critical en- 
tropy density (i.e. the entropy per lattice site) Sc, be- 
low which antiferromagnetic order emerges. Neverthe- 
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FIG. 1: (Color online) Temperature dependence of the en- 
tropy of the spin-1/2 quantum Heisenberg antiferromagnet 
and ferromagnet on the simple-cubic lattice. For both cases, 
results from quantum Monte Carlo simulations are shown for 
finite systems with L = 8 (dashed lines) and L — 16 (solid 
lines). Error bars for the shown data are below the line thick- 
ness. Vertical lines show the transition temperatures Tjv and 
Tc for the antiferromagnet and the ferromagnet, respectively. 
Green dashed-dotted curves denote fits to the algebraic low-T 
scaling S oc T 3 (S cx T 3 / 2 ) for the antiferromagnet (ferromag- 
net). 



less, thus far no systematic quantum Monte Carlo study 
concerning the entropy has been reported. On the other 
hand, different estimates of Sc for the Heisenberg anti- 
ferromagnet have been put forward recently. In Ref. [U, 
a Schwinger boson approach^* was employed, leading to 
the estimate that Sc is about 50% of the mean-field value 
S(j F = ks ln(2). A more recent study^ reported a value 
of Sc obtained from a fluctuation-corrected mean-field 
approach, leading to a reduction in Sc of only 20% from 
the mean- field value S C IF . This rather large discrepancy 
among current estimates of Sc calls for a clarification 
based on numerically exact simulations of the Heisen- 
berg model. Here, we present an unbiased estimate of 
the temperature dependence of the entropy and in partic- 
ular the critical entropy density by combining extended 
ensemble quantum Monte Carlo simulations with a finite- 
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size scaling analysis based on the known critical proper- 
ties of the quantum Heisenberg model. For comparison, 
we consider in the following both the ferromagnetic and 
the antiferromagnetic spin-1/2 Heisenberg model on the 
three-dimensional simple-cubic lattice. These models are 
described by a Hamiltonian 



H 



J ^ Si ■ Sj , 



(1) 



where Si denotes a spin-1/2 degree of freedom on lattice 
site i, and J is nearest- neighbor exchange interaction, 
with J > (J < 0), for the antiferromagnetic (ferromag- 
netic) case. From previous studies, the locations of the 
transition temperatures to the ordered phases have been 
determined as fcsTjv/J = 0.946(1)^ (Neel temperature) 
for the antiferromagnet, and k B T c /\J\ = 0.839(1)^ 
(Curie temperature) for the ferromagnet, based on quan- 
tum Monte Carlo simulations. 



II. METHOD 

In order to extract the temperature dependence of 
the entropy, we use an optimized extended ensemble 
approac h 17 ! 18 , that is based on a generalization of the 
Wang-Landau 19 algorithm to the case of quantum Monte 
Carlo simulation o 16 ' 20 , performed within the stochastic 
series expansion representational using multi-cluster de- 
terministic loop updates 22 . Within this approach, one 
obtains Monte Carlo estimates of the expansion coef- 
ficients g{n) of the high-temperature series expansion 
of the partition function Z in the inverse temperature 
p = l/(k B T), 



Z = Tre-? H = J29(n)l3 n , 



(2) 



n>0 



for a given system of N = L 3 lattice sites. Here, L de- 
notes the linear size of the simulated finite cube, and we 
employ periodic boundary conditions in all directions. 
From the expansion coefficients g(n), the free energy 



F=-lmz 



-iln£ 5 (n)/3", 



J!>0 



the internal energy 



ri>0 



and the entropy 



s = 



E-F 
T 



(3) 



(4) 



(5) 



are obtained as continuous functions of the temperature 
T. In practice, as discussed in Ref. [2(| we obtain the ex- 
pansion coefficients up to a upper cutoff A, that is cho- 
sen sufficiently large in order to reliably calculate the 



0.6 



0.55 



a 



0.5 



0.45 



0.4, 



, v ferromagnet 



. antiferromagnet 



0.05 



0.1 



0.15 
1/L 



0.2 



0.25 



FIG. 2: (Color online) Finite-size scaling of the critical en- 
tropy density Sc{L) with the inverse linear system size 1/L 
for the spin-1/2 quantum Heisenberg antiferromagnet and fer- 
romagnet on the simple-cubic lattice. The quantum Monte 
Carlo data are shown by symbols, and the dashed lines rep- 
resent fits of the finite-size data for L > 6 to the finite-size 
scaling form in Eq. (6). The results for Sc in the thermody- 
namic limit are indicated by symbols as well. 



thermodynamic properties of the finite system based on 
the coefficients g{n), n = 0, A down to a given tem- 
perature scale (3\- The required A scales linear in both 
/3\ and the system size N. In our simulations, we con- 
sider finite systems up to L = 20, for which A = 16000 
is required in order to obtain thermodynamic properties 
down to a temperature of about k B T « 0.7| J\. Repeated 
simulations with independent random streams results in 
statistical error bars on the thermodynamic quantities. 



III. RESULTS 

In Fig. 1, we present the results for the entropy from 
our simulations for both the antiferromagnetic and the 
ferromagnetic Heisenberg model. For both cases, we 
show data for systems with L = 8 and L = 16. Out- 
side a narrow range around the transition temperature, 
the entropy curves obtained for the two different system 
sizes compare well, and thus essentially represent the be- 
havior of the entropy density in the thermodynamic limit. 
Also shown in Fig. 1 are fits of the numerical data to the 
low-T scaling S oc T 3 and S oc T 3 / 2 for the antiferr omag- 
net and ferromagnet, respectively. This difference in the 
low-T scaling behavior relates to the linear (quadratic) 
dispersion of the magnetic spin-wave excitations in the 
antiferromagnetic (ferromagnetic) case^. 

However, in order to extract the critical entropy Sc in 
the thermodynamic limit, we need to perform a finite- 
size scaling analysis of the finite-size data Sc{L) taken 
at the transition temperature. For this purpose, we use 
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a finite-size scaling ansatz 

Sc-Sc(L) oci 1 ^, (6) 

where v is the critical exponent of the correlation length, 
and d the dimensionality of the system. Starting from 
S = —dF/dT, this relation follows from the scaling 
behavior of the free energy F oc with the corre- 
lation length £, that scales as £ oc \t\~ v when the re- 
duced temperature t — > 0, by the usual finite-size scal- 
ing ansatz, that at criticality £ relates to the linear sys- 
tem size L. For the current case of a three-dimensional 
system (d = 3) with 0(3) symmetry, v = 0.7112(5) 23 . 
Figure 2 shows that the finite-size data indeed fits well 
to the above scaling form in the range of system sizes 
L > 6. From the extrapolated values, we obtain 
S c = 0.341(5)fc B = 0.492(7)5$°" for the antiferromagn- 
tic model, whereas for the ferromagnet, a larger value of 
0.401(5)/c B = 0.578 (7) S% F results. We also performed 
an algebraic fit without biasing the exponent, which re- 
turned identical values of Sc within the statistical error 
bars. 

While the ferromagnet has a critical temperature Tq 
that is about 11% lower than T/v, it orders already below 
an entropy that is about 17% higher than the critical en- 
tropy of the antiferromagnet. The higher critical entropy 
directly exhibits a reduced efficiency of thermal disorder 
in destroying long-range magnetic order. The lower tran- 
sition temperature of the ferromagnetic is due to a sig- 
nificantly enhanced production of entropy upon heating 
up the system, which follows from both Fig. 1, as well as 
from the mentioned difference in the low-T behavior of 
the entropy. 



IV. CONCLUSIONS 



Based on numerically exact quantum Monte Carlo sim- 
ulations, we obtained the temperature dependence of the 
entropy for both the antiferromagnetic and the ferromag- 
netic spin-1/2 Heisenberg model on the simple-cubic lat- 
tice. Focusing on the region close to the magnetic or- 
dering transition, we employed a finite-size scaling anal- 
ysis to extract the critical entropy density below which 
long-ranged magnetic order sets in. In case of the anti- 
ferromagnetic model, our result for the critical entropy 
is consistent with a 50% reduction from the mean-field 
value, that was estimated recently from a Schwinger bo- 
son calculation^. This clarifies among the different values 
of the critical entropy reported previously^. In com- 
parison, we found the ferromagnet to be more robust to 
thermal disorder, which at criticality is about 17% larger 
for the ferromagnet than for the antiferromagnet in terms 
of the critical entropy. We consider the current results 
for the temperature dependence of the magnetic entropy 
of value also for future calculations of the entropy for 
the full fermionic Hubbard model beyond the strongly 
interacting regime. 
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